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Abstract 

A simple model based on relativistic geometry and final-state hadronic rescattering is used to 
predict pion source parameters extracted in two-pion femtoscopy studies of proton-proton collisions 
at ^/s = 200 GeV. From studying the momentum and particle multiplicity dependences of these 
parameters in the context of this model and assuming a very short hadronization time, flow-like 
behavior is seen which resembles the flow behavior commonly observed in relativistic heavy-ion 
collisions. 
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In the field of fluid dynamics, where a fluid can be either a liquid or a gas, fluid flow 
is generally defined as the transport of a certain amount of fluid mass or volume across 
a surface in a given time interval. This concept of fluid flow has been adopted by the 
relativistic heavy-ion collision community to describe the behaviors of some observables 
seen in experiments which appear to be flow-like in nature . The justification for 

having a fluid dynamics picture of a relativistic heavy-ion collision is that in these collisions, 
for example colliding beams of Au + Au at 200 GeV per nucleon pair in the center-of-mass 
frame, thousands of particles such as partons, e.g. quarks and gluons, hadrons, e.g. nucleons 
and pions, and leptons, e.g. electrons and muons, participate and are generated out of the 
vacuum in a violent collision. Thus the concept of characterizing the "bulk properties" of 
such collisions, for example the dynamics of the size and shape of the interaction region, by 
a relativistic fluid flow where the fluid is made up of thousands of particles seems reasonable 
to some approximation. On the other hand, if one were to try to impose a fluid dynamics 
picture on proton-proton collisions at similar energies, e.g. p + p collisions at ^/s = 200 
GeV, in which typically much fewer than 100 particles participate and are produced in the 
collision, such an approximation would seem unreasonable due to the paucity of particles 
available to compose the "fluid." 

The goal of the present work is to use a simple model to study one of these bulk prop- 
erties in p + p collisions at ^/s = 200 GeV, namely the dynamics o: 
of the interaction region, using the method of two-pion femtoscopy 
Hanbury-Brown-Twiss interferometry (or HBT, which has also been used to measure the 
sizes of stars . The HBT observables, which are radius parameters, have been shown 
in Au + Au collisions at ^/sN^ = 200 GeV to exhibit radial flow-like behavior in their de- 
pendences on the pion momenta and the particle multiplicity of the collision, such that the 
radius parameters decrease for increasing pion momenta (higher velocity fluid with smaller 
correlation length) and increase for increasing particle multiplicity (higher density fluid pro- 
ducing more expansion) [8] . If the model predicts that there are similar dependences also 
present in p + p collisions at y/s = 200 GeV it would help to establish the possibility of 
flow-like effects in these collisions as well as characterize the mechanism by which these 
effects are generated. 

The model calculations are carried out in four main steps: 1) simulate p+p collisions using 
the standard collision-generator code PYTHIAjol, 2) employ a simple space-time geometry 



the size and shape 
^, B] also known as 
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picture for the hadronization of the PYTHIA-generated hadrons, 3) calculate the effects of 
final-state rescattering among the hadrons, and 4) calculate the HBT radius parameters. 
These steps will now be discussed in more detail. 

n 

The p + p collisions were modeled with the PYTHIA code, version 6.409 [9]. The parton 
distribution functions used were the same as used in Ref. [llj. Events were generated in 
"minimum bias" mode, i.e. setting the low-pT cutoff for parton-parton collisions to zero and 
excluding elastic and diffractive collisions. To obtain good statistics for this study 5 x 10^ 
events were generated with ^/s = 200 GeV. Information saved from a PYTHIA run for 
use in the next step of the procedure were the momenta and identities of the "direct" (i.e. 
redundancies removed) hadrons (all charge states) vr, K, p, n, A, p, u, r], r]', 0, and K*. 
These particles were chosen since they are the most common hadrons produced and thus 
should have the greatest effect on the hadronic observables in these calculations. 

The space-time geometry picture for hadronization from a. p + p collision in the model 
consists of the emission of a PYTHIA particle from a thin uniform disk of radius 1 fm 
in the x — y plane followed by its hadronization which occurs in the proper time of the 
particle, r. The space-time coordinates at hadronization in the lab frame {xh,yh, Zh,th) for 
a particle with momentum coordinates {Px,Py,Pz), energy E, rest mass mo, and transverse 
disk coordinates (xo,yo), which are chosen randomly on the disk, can then be written as 



Xh = Xo + T— (1) 

mo 

yH = yo + (2) 
mo 

Zh = r— (3) 

mo 

= r— (4) 

mo 

The simplicity of this geometric picture is now clear: it is just an expression of causality 
with the assumption that all particles hadronize with the same pro per time, r. A sim- 
ilar hadronization picture has been applied to e+ — collisions 
collisions [11]. For all results presented in this work, r will be set to 0.1 fm/c to be con- 
sistent with the results found in the Tevatron study which had reasonable agreement with 
measurements. 

The hadronic rescattering calculational method used is similar to that employed in previ- 
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ous studies js], [ll| . Rescattering is simulated with a semi-classical Monte Carlo calculation 
which assumes strong binary collisions between hadrons. Relativistic kinematics is used 
throughout. The hadrons considered in the calculation are the most common ones: pions, 
kaons, nucleons and lambdas (vr, K, N, and A), and the p, uj, rj, rj', (p, A, and K* resonances. 
Although A formation was included in the rescattering process due to its large formation 
cross section in meson-baryon and baryon-baryon rescattering, they were not input directly 
from PYTHIA since they decay quickly and relatively few are initially produced in these 
collisions. For simplicity, the calculation is isospin averaged (e.g. no distinction is made 
among a vr^, tt^, and tt^). Starting from the initial stage (t = fm/c), the positions of all 
particles in each event are allowed to evolve in time in small time steps (At = 0.1 fm/c) 
according to their initial momenta. At each time step each particle is checked to see a) 
if it has hadronized {t > th, where th is given in Eq. (4)), b) if it decays, and c) if it is 
sufficiently close to another particle to scatter with it. Isospin-averaged s-wave and p-wave 



cross sections for meson scattering are obtained from Prakash et al. 12| and other cross sec- 

n 

tions are estimated from fits to hadron scattering data in the Review of Particle Physics[13l|. 
Both elastic and inelastic collisions are included. The calculation is carried out to 50 fm/c 
which allows ample time for the rescattering to finish. Note that when this cutoff time is 
reached, all un-decayed resonances are allowed to decay with their natural lifetimes and their 
projected decay positions and times are recorded. The validity of the numerical methods 
used in the rescattering code have recently been studied using the subdivision method, the 
results of which have verified that the methods used are valid 

In carrying out a two-pion HBT study the two-pion coincident count rate is calculated 
along with the one-pion count rates for reference. The two-pion correlation function for pions 
binned in momenta pi and p2, C*(pi, P2), is constructed event-by-event from the coincident 
countrate, iV2(pi,p2) and one-pion countrate summed over events, A^i(p), as 

ryf \ ^2(pi,P2) 

A^i(pi)A^i(P2) 

It is usually convenient to express the six- dimensional C(pi, P2) in terms of the four- vector 
momentum difference, Q = |pi — P2I by summing Eq.([5]) over momentum difference, 

r^(n\ - ^2(Pi,P2) _ A{Q) 

P1,P2(Q) 

where A{Q) represents the "real" coincident two-pion countrate and B{Q) the "background" 



two-pion countrate composed of products of the one-pion countrates, all expressed in Q. 
In practice, B{Q) is the mixed event distribution, which is computed by forming pion-pairs 
from different events. 

For the HBT calculations from the model, a three-dimesional two-pion correlation func- 
tion is formed using Eq.([2]) and a Gaussian function in momentum difference variables is 
fitted to it to extract the pion source parameters similar to what is done in experiments 
Boson statistics are introduced after the rescattering has finished using the standard method 
of pair-wise symmetrization of bosons in a plane- wave approximation The three- 

dimensional correlation function, C {Q side, Q out, Q long), is then calculated in terms of the 
momentum-difference variables Qside, which points in the direction of the sum of the two 
pion momenta in the transverse plane, Qout, which points perpendicular to Qside in the 
transverse plane and the longitudinal variable along the beam direction Qiong- 

The final step in the calculation is extracting fit parameters by fitting a Gaussian param- 
eterization to the model-generated two-pion correlation function given by, Gj 



C(yQside, Qout, Qlong) (7) 

/ long-^longl 



G[l + Aexp( Qlide^Mde Qout^out Q'lonq^'lon 



where the i?-parameters, called the radius parameters, are associated with each momentum- 
difference variable direction, G is a normalization constant, and A is the usual empirical 
parameter added to help in the fitting of Eq. ([7]) to the actual correlation function (A = 1 
in the ideal case). The fit is carried out in the conventional LCMS frame (longitudinally 
comoving system) in which the longitudinal pion pair momentum vanishes 6|]. 

Figure [T] shows a sample three-dimensional correlation function from the model projected 
onto the Qside axis with projected fit to Eq. ([71). The other variables, Qout and Qiong, are 
integrated up to 0.02 GeV/c. The cuts on the pion momenta are —0.5 < y < 0.5, where y is 
rapidity, 0.15 < px < 0.8 GeV/c, where px is transverse momentum, and 0.15 < < 0.25 
GeV/c (see below). These cuts in y and pr are used throughout to duplicate those used in 
experiments. Although only a small fraction of the full correlation function is shown in this 
plot, it gives some idea of the quality of the Gaussian fit to the model. For Qside < 0.05 
GeV/c the fit is seen to undershoot the model which is a common feature also seen in 



experiments ioi p + p collisions ly, [iTj]. The A parameter for this fit is 0.330, which is far 



less than the ideal value of unity, another feature commonly seen in p + p experiments. The 
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FIG. 1: Sample two-pion correlation function with three-dimensional Gaussian fit (Eq.Q) pro- 
jected onto the Qside axis from the Model. The other variables, Qout and Qiong, are integrated up 
to 0.02 GeV/c. The cuts on the pion momenta are —0.5 < y < 0.5, 0.15 < px < 0.8 GeV/c, and 
0.15 <kT < 0.25 GeV/c. 



source of this non-Gaussian and non-ideal behavior in the model is the presence of medium 
to long-lived resonances such as the u, t], and r^' which give a component to the correlation 
function representing a much larger pion source than the majority of the pions in the collision 
and thus producing the narrower shape in the correlation function near Q = 0. 

Figure [2] shows the pion momenta and particle multiplicity dependences of the radius 
parameters from the model, both for the full model calculations (top plots) and for model 
calculations in which hadronic rescattering is turned off (bottom plots). The pion momenta 
and particle multiplicity are represented by the quantities kx = (pti + Pt2)/2 which is 
the average transverse pion momentum of the pair, and {1/N)dn/ drj evaluated at = , 
which is the (pseudo)rapidity density of charged particles at mid-rapidity, respectively. The 
dashed lines are fits to the model points, linear for the fcy-dependence plots and logarithmic 
for the (l/A^)(in/(i?7-dependence plots. Focusing on the top two plots first, the full model 
calculations, it is clear that the radius parameters have the same qualitative dependences on 
kx and rapidity density as observed in heavy-ion collisions, namely decreasing with increasing 
krp and increasing with increasing rapidity density (i.e. particle multiplicity). Another trend 
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FIG. 2: kT and rapidity density dependences of the radius parameters from the model. The top 
plots show the full model calculations and the bottom plots show model calculations in which 
hadronic rescattering is turned off. The dashed lines are linear and logarithmic fits to the model 
points for the /ct and rapidity density plots, respectively. 



seen in Figure [2] which is also observed in heavy-ion collisions is Riong > Rout > Rside- That 
hadronic rescattering is the main source of these effects in the model is seen by comparing the 
top plots with the bottom plots for which rescattering is turned off in the model. In addition 
to the overall scales of the radius parameters being significantly smaller with rescattering 
turned off, it is also seen that all of the dependences which were seen in the top plots are 
either greatly reduced or absent in the bottom plots. A small degree of the kx dependence 
is seen to linger for R^ut and Riong which is mainly caused by the resonances present in the 
model. 

In order to make a more quantitative comparison between the dependence seen in 
Figure [2] in the full model for p + p collisions and that measured in heavy-ion collisions, a 
calculation has been made of the ratios of the radius parameters extracted in the STAR 
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experiment for Au + Au collisions at ^/snn = 200 GeV[8| to those in Figure [2l The same 
kinematic cuts on the pions as used in Figure [2] were used by the STAR experiment and the 
STAR results were extracted for a centrality cut of 0-5%. Figure [3] shows a plot of these 
ratios vs. kx- The STAR radius parameters used in calculating the ratios are also plotted. 
Comparing the kj- dependence of the model p + p in the upper left plot of Figure [2] with the 
STAR radius parameters shown in Figure [3] the qualitative similarity of the plots is evident. 
As seen, the ratios are approximately flat in kx, but with a slight increasing tendency hinting 
that the decrease of the p + p radius parameters with kx from the model is slightly stronger 
than that measured in the Au + Au data. Preliminary measurements have been made of 
radius parameters from p + p colhsions at y/s = 200 GeV by STAR and experimental ratios 
have been calculated as in Figure [3] showing a similarly fiat dependence on kr [3] • 

It is difficult to make a quantitative comparison between the model p + p results for the 
rapidity density shown in Figure [2] and heavy-ion experiments in a similar way as was done 
for the kr dependence since the range in rapidity density needed to represent the Au + Au 
measurements is about 100-650, which is much larger than that seen in Figure [2] of 0.3-10. 
A qualitative comparison can be made in that it can easily be shown that the rapidity 
density dependence of STAR Au + Au radius parameters from ^/sN^ = 200 GeV collisions 
is approximately logarithmically increasing with increasing rapidity density as is seen in 
Figure [10, 19]. 

In conclusion, it has been shown that flow-like effects observed in relativistic heavy-ion 
collision experiments can be reproduced in p+p collisions at y/s = 200 GeV by a simple model 
based on relativistic geometry and final-state hadronic rescattering with a short proper time 
for hadronization of 0.1 fm/c. In the model, the fiow-like effects are driven by the hadronic 
rescattering which in turn is sensitively controlled by the hadronization time which sets the 
initial particle density. If a long hadronization time were used, e.g. r = 1 fm/c, the initial 
particle density would be low and little rescattering would take place such that the full 
model results would more resemble the bottom plots in Figure El Thus to the extent of the 
agreement shown above between the present model calculations and experiment, another 
implicit result of this study is that the hadronization time in p + p collisions at this energy 
appears to be very short. 

The author wishes to acknowledge financial support from the U.S. National Science Foun- 
dation under grant PHY-0653432, and to acknowledge computing support from the Ohio 
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FIG. 3: Ratios of STAR Au + Au radius parameters to model p + p radius parameters vs. kj-- The 
STAR radius parameters used in calculating the ratios are also plotted. The horizontal solid lines 
are the averages over of the ratios for a given radius parameter. The dashed lines are linear fits 
to the STAR radius parameters. 



Supercomputing Center. 



|humanic@mps.ohio-state.edu| 

I. Arsene et al. [BRAHMS collaboration], Nucl. Phys. A 757, 1 (2005). 
B. B. Back et al. [PHOBOS collaboration], Nucl. Phys. A 757, 28 (2005). 
J. Adams et al. [STAR collaboration], Nucl. Phys. A 757, 102 (2005). 
K. Adcox et al. [PHENIX collaboration], Nucl. Phys. A 757, 184 (2005). 
T. J. Humanic, Int. J. Mod. Phys. E 15, 197 (2006). 

M. A. Lisa, S. Pratt, R. Soltz and U. Wiedemann, Ann. Rev. Nucl. Part. Sci. 55, 357 (2005). 

R. Hanbury Brown and R. Q. Twiss, Nature 177, 27 (1956). 

J. Adams et al. [STAR Cohaboration], Phys. Rev. C 71, 044906 (2005). 



T. Sjostrand, L. Lonnblad, S. Mrenna and P. Skands, arXiv:hep-ph/0603175 (March 2006). 
T. Csorgo and J. Zimanyi, Nucl. Phys. A 512, 588 (1990). 
T. J. Humanic, Phys. Rev. C 76, 025205 (2007). 



9 



[12] M. Prakash, M. Prakash, R. Venugopalan and G. Welke, Phys. Kept. 227, 321 (1993). 

[13] W. M. Yao et al. [Particle Data Group], J. Phys. G 33, 1 (2006). 

[14] T. J. Humanic, Phys. Rev. C 73, 054902 (2006). 

[15] T. J. Humanic, Phys. Rev. C 34, 191 (1986). 

[16] T. Alexopoulos et al., Phys. Rev. D 48, 1931 (1993). 

[17] J. L. Bailly et al. [NA23 Collaboration and EHS-RCBC Collaboration], Z. Phys. C 43, 
(1989). 

[18] Z. Chajecki [STAR Collaboration], Nucl. Phys. A 774, 599 (2006). 

[19] B. B. Back et al, Phys. Rev. Lett. 91, 052303 (2003). 



10 



